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Four-point functions and the permutation group 5*4 


Gernot Eichmann, Christian S. Fischer, and Walter Heupel 
Institut fiir Theoretische Physik, Justus-Liebig-Universitat Giessen, 35392 Giessen, Germany. 

Four-point functions are at the heart of many interesting physical processes. A prime example is 
the light-by-light scattering amplitude, which plays an important role in the calculation of hadronic 
contributions to the anomalous magnetic moment of the muon. In the calculation of such quantities 
one faces the challenge of finding a suitable and well-behaved basis of tensor structures in coordinate 
and/or momentum space. Provided all (or many) of the external legs represent similar particle 
content, a powerful tool to construct and organize such bases is the permutation group S 4 . We 
introduce an efficient notation for dealing with the irreducible multiplets of Si, and we highlight the 
merits of this treatment by exemplifying four-point functions with gauge-boson legs such as the four- 
gluon vertex and the light-by-light scattering amplitude. The multiplet analysis is also useful for 
isolating the important kinematic regions and the dynamical singularity content of such amplitudes. 
Our analysis serves as a basis for future efficient calculations of these and similar objects. 


I. INTRODUCTION 

In the study of the structure of elementary and com¬ 
posite particles one often faces the situation that the 
physical properties of interest are encoded in higher 
n—point functions. These appear in different contexts. 
For example, hadron properties are experimentally ex¬ 
tracted from nucleon-lepton or NN scattering, tttt scat¬ 
tering, Compton scattering or pion electroproduction 
amplitudes, which are all four-point functions. The struc¬ 
ture of electrons is inferred from Compton scattering and 
even the structure of photons can be measured in prin¬ 
ciple via light-by-light (LbL) scattering on the photon 
mass shells. The corresponding off-shell photon four- 
point function, although not directly measurable, is a 
particularly interesting case since it receives hadronic cor¬ 
rections that play an important role in the muon {g — 2) 
puzzle [1]. On a technical level, the photon four-point 
function is challenging because of the gauge boson nature 
of its external legs, but it is also subject to potentially 
extensive simplifications due to the corresponding Bose 
symmetry. 

Higher n—point functions are also important at the 
level of quarks and gluons: hadrons appear as poles in 
the qq four-point and qqq six-point correlation functions, 
whose residues define the qq and qqq Bethe-Salpeter 
vertex functions, and their properties are linked with 
the fundamental quark-gluon, three-gluon and four-gluon 
vertices. With respect to the recent experimental discov¬ 
eries of a variety of potential exotic states (XYZ-states) 
in the heavy quark region, a further interesting quantity 
is the qqqq eight-point function, with tetraquark ampli¬ 
tudes as residues on the corresponding pole positions. 
Tetraquark amplitudes are five-point functions with four 
(anti-)quark legs and display symmetries similar to the 
corresponding four-point functions. 

It is desirable to develop the technology to deal with 
such quantities efficiently and make the most possible 
use of the underlying symmetries. A generic n—point 
function depends on n — 1 independent momenta and 


has the form 

N 

=Y^fi{...)Tr'"{PU---Pu) ( 1 ) 

i 

in momentum space, modulo potential flavor and color 
factors, where /i, v denote Lorentz and/or Dirac indices. 
It can be decomposed into N Dirac-Lorentz tensors ' 
with corresponding dressing functions fi. The fi depend 
on M = n(n — I)/2 — m independent Lorentz invariants 
Tij P 21 Pi ' P 2 ■ ■ ■ that constitute the phase space, where 
m is the number of legs that are on-shell. For example, 
in linear covariant gauges one finds 

• (TV, M) = (14, 3) for the three-gluon vertex [2], 

• {N,M) = (136,6) for the four-gluon vertex and, 
without implementing gauge invariance, also for the 
LbL amplitude (see Sec. VIA), 

• (TV, M) = (64, 5) for a nucleon Bethe-Salpeter am¬ 
plitude [3] and (128,5) for that of a A—baryon [4], 

• and (TV, M) = (256,9) for a tetraquark Bethe- 
Salpeter amplitude [5]. 

Clearly, for a growing number of external legs this be¬ 
comes unmanageable and poses a challenge for theoret¬ 
ical approaches. For example, in the non-perturbative 
region of QCD the fundamental two-point functions are 
reasonably well understood by now, partially owing to 
their relatively simple structure. Progress in determin¬ 
ing three-point functions is underway, see e.g. [6-9] and 
references therein. Concerning four- and higher n—point 
functions, however, only exploratory studies exist so far 
[10-12]. Within QED, an important and prominent ex¬ 
ample is the photon four-point function, which received 
an increasing amount of interest recently in connection 
with the problem of the anomalous magnetic moment of 
the muon. Its hadronic corrections account for a sizable 
contribution to the total error budget of the theoretical 
determination of {g — 2)^. With efforts for a substantial 
decrease of the experimental error underway [13, 14] var¬ 
ious theoretical approaches to capture the rich physics 
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of the photon four-point function on a quantitative level 
have been explored [15-22]. Clearly, a suitable and well- 
organized basis of tensor structures of this object will 
greatly facilitate these studies. 

In general, the organizing principles for a basis of ten¬ 
sor structures of an n-point function are gauge invariance 
(if gauge-boson legs are involved), momentum counting 
and permutation-group symmetries. Correlation func¬ 
tions that are subject to gauge-invariance constraints can 
be split into a ‘gauge’ part and a transverse contribution, 
where the gauge part is often used as a zeroth-order ap¬ 
proximation to the full result. For gange-invariant quan¬ 
tities the gauge part vanishes, which complicates matters 
further because the transverse part is also subject to an- 
alyticity constraints. The momentum counting entails 
that tensor structures with higher momentum powers are 
less important, and in practice it is often sufficient to 
restrict oneself to the simplest momentum-independent 
tensors to obtain a reasonable approximation of the full 
n-point function. Finally, permutation-group symmetries 
are useful because they allow one to arrange both tensor 
structures and momentum variables into multiplets and 
thereby isolate the most relevant momentum dependence 
of the dressing functions fi. 

While our main motivation for the present paper is 
the muon g — 2 problem, we will not repeat the discus¬ 
sion of the underlying physics problems here. Instead, 
we will focus on technical aspects: we develop the termi¬ 
nology to deal with permutation-group symmetries, and 
we apply it to investigate the phase space in four-point 
amplitudes and construct appropriate tensor bases for 
them. To some extent our results are also applicable to 
other systems where permutation-group symmetries play 
a role, such as the four-gluon vertex or the tetraquark 
Bethe-Salpeter amplitude, and in principle they can be 
generalized to higher n—point functions. 

The paper is organized as follows. We exemplify the 
basic problems in constructing gauge-invariant tensor 
bases in Sec. II for a scalar two-photon current. In Sec. Ill 
we illustrate our multiplet notation for the permutation 
group S 3 and its application to baryon flavor wave func¬ 
tions. We subsequently generalize it to the permutation 
group S 4 in Sec. IV, where we discuss the color struc¬ 
ture of the four-gluon vertex as an example. In Sec. V 
we apply the terminology to examine the phase space in 
four-point functions, and we finally construct appropri¬ 
ate tensor bases for vector four-point functions in Sec. VI. 
We summarize and conclude in Sec. VII. 


II. EXAMPLE: SCALAR TWO-PHOTON 
CURRENT 

Ultimately our goal is to establish tensor bases for vec¬ 
tor four-point functions that implement the constraints 
from Bose symmetry, gauge invariance and analyticity. 
The template for this is the analogous case of Compton 
scattering that has been investigated in Refs. [23-29]. We 
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FIG. 1. Kinematics for a scalar two-photon current. 


will discuss four-point functions in Sec. VI, but to illus¬ 
trate the main points we first consider a simpler system, 
namely a two-photon current with scalar (0"'"“'") quantum 
numbers. 

The current r^'^(fc, Q) depends on two independent 
momenta (see Fig. 1); we denote the outgoing and in¬ 
coming photon momenta by 

ki=k+^, k2 = k-^, (2) 

where Q = ki — k 2 is the total bound-state momentum 
and k = (ki + k 2)/2 is the average photon momentum. 
Electromagnetic gauge invariance entails that r^'^(fc, Q) 
is transverse with respect to and : 

kV YV’' = 0, F^”^ = 0 . (3) 

Since the current has two photon legs, transversality and 
analyticity require it to be at least quadratic in the pho¬ 
ton momenta. 

The most general decomposition of the current consists 
of five tensor structures, which we can express in terms 
of the photon momenta fci and k 2 '. 

5 

= E = /i +/2 fcr 

i=i + /3 e, k’i + u {k^i fcr + k^ k-2 ) 

+ f5^{ki fci - k^) ■ 


The dressing functions or form factors fi{k‘^,uj,Q‘^) de¬ 
pend on the Lorentz invariants /c^, and tv = k ■ Q, ot 
equivalently 


= k'^ + ^ ±k-Q = r]+±u}, 
ki ■ k2 = k"^ - ^ = r]-. 


( 5 ) 


Bose symmetry entails T^''{k,Q) = T''^{—k,Q). The 
Lorentz invariants k^ and are symmetric whereas w 
is antisymmetric. To ensure Bose symmetry at the level 
of the basis elements, the tensor structure for /s includes 
an angular prefactor w so that all form factors fi are even 
in w, so they depend on k^, and In other words, 
all basis elements and all form factors are now singlets 
under the permutation group 82 - The consequences of 
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analyticity are not yet manifest in this basis; the require¬ 
ment we will need is that the fi are nonsingular in any 
kinematic limit. 

Instead of working out Eqs. (3) it is more convenient in 
practice to equate the transverse projection of the current 
with the current itself: 

rjifia pa/3 rpfSi' ^ p/ij^ ^ 0 ^ 

Here, = 6 ^'' — /cf k'^/kf is a transverse projector 
with respect to photon i. Only the tensor structures for 
fi and /2 in Eq. (4) survive the projection (we will call 
them transverse survivors^ in the following) whereas the 
others are longitudinal. This leads to three conditions 
for the fi whose solution is 

fi = -V-f2 + (??+ 

f3 = V-f5, (7) 

/4 = -V+f5 , 


with ri± defined in Eqs. (5). Here we exploited the con¬ 
straint that the fi should be nonsingular: had we solved 
for any other set of three form factors, we would have 
introduced denominators that can become singular, con¬ 
trary to our initial assumption. Note also that the Bose 
symmetry of the fi was crucial: without the antisym¬ 
metric variable oj in the last row of Eq. (4) the resulting 
factors of oj would obscure the correct solution. Hence, to 
solve the transversality constraints we should start from 
a basis made of permutation-group singlets. 

When substituting the solution into Eq. (4) , the result 
for the transverse current becomes 

r- = / 2 h-,_n]- 

+ h [t-s + (?7+ - n + Ji-T-s - V+ta] , 

These two new tensor structures do not look particularly 
enlightening, but they can be cast in a compact form if 
we define 


= k, ■ kj ,5^" - fc; kf . (9) 

This quantity is transverse with respect to k^ and klf 
without introducing a pole, and for i = j li becomes 
proportional to the usual transverse projector: = 

kfTf' with Tf' defined below Eq. ( 6 ). The resulting 
current in Eq. ( 8 ) can now be written as 

=-f 2 t^l 2 +ht'tlt^ 2 , ( 10 ) 

from where the transversality and analyticity properties 
can be read off directly.^ The two transverse structures 


^ Following Tarrach’s procedure for Compton scattering [24], we 
would write instead of Eq. (6): 


rj-iftCX 


~ h - . h ■ ' 

rb j 


( 11 ) 


which leads to the same result. 


exhibit the correct power counting. Since the current can 
be expressed entirely in terms of photon momenta, the 
photon momentum powers are equivalent to their mass 
dimension: has mass dimension two and tftf has 

dimension four. The transverse tensors read explicitly: 

= ki-k2 - klf k'i, 

C t22 = - k^, kZ) {kl - k^ k!f) (^2) 

— h’' 

— ^2 ^ ^2 A /^2 ^2 

+ ki-k 2 kf k^. 


Their form factors /2 and /a are free of kinematic singu¬ 
larities and zeros and become constant in any kinematic 
limit of vanishing photon momenta. Hence, their only 
possible singularities are of dynamical origin. The ba¬ 
sis ( 10 ) is therefore ‘minimal’. 

This simple example highlights the main features that 
we will again encounter in the discussion of the LbL am¬ 
plitude. What if the two-photon current happened to 
be not transverse, for example as a consequence of an 
incomplete dynamical calculation that violates gauge in¬ 
variance? We can add again the tensor structures for 
those coefficients in Eq. (7) that we found to be linearly 
dependent: 

= g, <5^-^ + g, k>^ klf + 34 (kf fcr + k!f k!f ) 

f 4^"^ I f 4-ai' 

/2 ti2 + 75 Cll ^22 ■ 

Violations of gauge invariance will then become manifest 
in nonzero functions gi, 33 and 34 . If we apply another 
transverse projection to Eq. (13) and express the result in 
the ‘physical’ basis t^lf and ^22 ) these gauge artifacts 
will produce unphysical singularities in the two-photon 
current. 


rji(-LOL j^Q(/3 _ y* 


h + 


91 

^2 1.2 
^X ^2 


,/j,a ,otv 

^Xl ^22 5 


(14) 


and therefore invalidate the extraction of form factors in 
the limits where the photon momenta vanish. Ultimately, 
in a fully gauge-invariant calculation all 3 / must cancel, 
so instead of projecting transversely we could simply omit 
the gi terms in Eq. (13). The form factors / 2 , /s of the 
transverse remainder would then still provide a sensible 
estimate of the result. 

Eq. (13) is the general decomposition of a scalar- 
vector-vector amplitude in terms of a ’gauge part’ and 
a transverse part that implements the analyticity con¬ 
straints. It also applies to a scalar glueball amplitude 
[30, 31], and analogous constructions are well-known for 
the quark-photon and quark-gluon vertex, the three- 
gluon vertex etc. In those cases, the gauge parts are 
nonzero and constrained by Ward-Takahashi or Slavnov- 
Taylor identities, and they are not artifacts but usually 
even the dominant contributions. For example, the gauge 
part in the quark-photon vertex is the Ball-Chiu ver¬ 
tex [32] that is known to give a large contribution to 
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hadronic form factors. Only for gauge-invariant ampli¬ 
tudes the gauge parts are forced to be zero, and in combi¬ 
nation with the power counting this leads to the problem 
of unphysical singularities if gauge invariance is not re¬ 
spected. 

The example also provides us with another, more direct 
way to arrive at Eq. (13). Since the transverse projection 
of Eq. (4) leaves only two survivors, it is clear that there 
can be only two transverse tensor structures. Using the 
definition (9), we could have immediately written down 
the two lowest-dimensional transverse tensors that are 
permutation-group singlets and free of kinematic singu¬ 
larities: t ^2 ^11 ^22 ■ When expressed in terms of the 

generic basis via Eq. ( 8 ), the former depends on T 2 and 
the latter on These are the only dependencies that 
do not involve any kinematic prefactors. Therefore, the 
only kinematically safe option is to eliminate the contri¬ 
butions from T 2 and in Eq. (4) in favor of the new 
transverse tensors, thereby leaving ti, and T 4 as the 
gauge part. The result is identical to Eq. (13). 

To facilitate the following discussion, we will refer to 
the two types of bases in Eq. (4) and Eq. (13) as ’type I’ 
and ’type IF, respectively: 

F = Fsur + F L = Fgauge + Fy . (15) 

I type II 

The type-I basis consists of a ‘survivor’ and a longitudinal 
part and the type-II basis of a gauge and a transverse 
part. The transverse survivors depend upon the type 
of projection: Eqs. ( 6 ) and (11) will produce different 
survivors, but their total number is the same. Although 
the gauge parts are usually called ‘longitudinal’, neither 
Fsur nor Fgauge are fully transverse or longitudinal. A 
transverse or longitudinal projection will rather result in 
the following structure: 

r(F)=T(F,ur)=T(Fgauge)+rT, 

L(F) = L(F,ur) + L{Tl) = T(rgauge) • ^ ’ 

Whereas the basis dimensions are identical, 

dim(Fsur) = dim(rT), dim(Fgauge) = dim(FL), (17) 

there is no one-to-one relation between their dressing 
functions, which is also apparent in the two-photon cur¬ 
rent example of Eqs. (4) and (13). 

The example also illustrates another generic problem: 
type-I bases are straightforward to write down, whereas 
the construction of type-II bases, in which the underly¬ 
ing physics becomes more transparent, is much more dif¬ 
ficult. QCD vertices that are contracted with transverse 
gluon propagators in Landau gauge have the form T(r) 
with Fgauge 7 ^ 0, hence type-I bases are usually sufficient 
for their analyses.^ In the case of gauge-invariant photon 


^ In the case of the quark-photon vertex, Eq. (75) in Ref. [33] con- 


amplitudes only Ft survives, which makes it mandatory 
to find an appropriate basis of type II. 

The program in applying these findings to the photon 
or gluon four-point functions is then the following: 

• Work out the irreducible representations of the per¬ 
mutation group S 4 . 

• Cast the Lorentz-invariant momentum variables 
into S '4 multiplets and work out the structure of 
the phase space, the relevant kinematic limits, etc. 

• For the photon four-point function, write down a 
generic Lorentz tensor basis of type I and cast it 
into S4 singlets. For the four-gluon vertex, ar¬ 
range also the color structures into «S '4 multiplets 
and combine them together with the Lorentz ten¬ 
sors into S4 singlets. 

• Apply the transversality constraints and derive a 
type-II basis which is made of S4 singlets. It should 
be the sum of a gauge and a transverse part, where 
the gauge part must vanish if the amplitude is 
gauge invariant. 

Since the construction of permutation-group multiplets 
requires a certain amount of formalism, we will first dis¬ 
cuss the simpler case of the permutation group S 3 and 
postpone its application to the group S4 to Sec. IV. 


III. PERMUTATION GROUP S3 
A. Multiplets 

The permutation group S3 consists of 3! = 6 elements. 
The group manifold is visualized by the Cayley graph in 
Fig. 2: any permutation of an object /123 can be recon¬ 
structed from a transposition P 12 and a cyclic permuta¬ 
tion P 123 . The former interchanges the indices I fA 2 and 
the latter is a cyclic permutation 1— 7 ^ 2 , 2—3—>■!. 
The group elements acting on /123 are given by 

1 , P12, 

P13 P12 = P123 , P23 = P12 P123 , ( 18 ) 

P23 P12 = P123 I P13 = P12 P123 

and can be visualized by paths along the Cayley graph, 
starting from the lower left corner. 


stitutes the type-I basis and Eq. (89), together with the Ball-Chiu 
vertex of Eq. (72), the type-II basis. The same applies to the 
quark-gluon vertex except that Fgauge follows from a Slavnov- 
Taylor identity that contains the quark-ghost kernel. For the 
three-gluon vertex. Table V in Ref. [7] is of type I and Eq. (66) 
is the type-II basis. 
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FIG. 2. Cayley graph for the permutation group S 3 . Any 
permutation can be reconstructed from a transposition P 12 
and a cycle P123. 


with the two-dimensional orthogonal representation ma¬ 
trices (M^ denotes the matrix transpose) 

Mi 2= ( J J), Mi23=^_j3 ^ )■ 

S is totally symmetric and invariant under any permuta¬ 
tion. A is totally antisymmetric under exchange of any 
two entries and therefore it picks up a minus sign under a 
transposition. Both subspaces are one-dimensional. The 
doublets form a two-dimensional subspace; they trans¬ 
form under the same matrix representations M 12 and 
Mi 23 from where all other ones can be reconstructed via 
Eq. (18), e.g.: P 23 —(M 12 Mi 23 )^. Their upper (lower) 
components are antisymmetric (symmetric) under trans¬ 
positions P 12 and we denote them by 


The six permutations of /123 can be rearranged into a 
symmetric singlet S, an antisymmetric singlet (‘antisin¬ 
glet’) A, and two doublets Vj {j = 1 , 2 ) whose explicit 
form can be chosen as^ 

S = ijt +^t. 

A = Tp- + ij)- A ip- , 


'^1 = 



(23) 


In practical applications the tensor products of S 3 
multiplets are often of interest. Given two sets of sin¬ 
glets 5, S', antisinglets A, A' and doublets V = [a, s], 
V = [o', s'], there are 16 possible combinations 

{S, A, a, s} X {S', A', a', s'} (24) 


Vi 

V 2 


V’2 -V’3 

^ {i^2 + ^3 - 2^r)' 


where 


(19) 


^f=P±fl23, i’t = P±f231. i’t=P±f312 ( 20 ) 

and P± = 1 ± P 12 . They transform under the irreducible 
representations of S 3 , which correspond to the following 
Young diagrams: 


which we can again arrange into multiplets. Of course 
the products of two singlets {SS') or antisinglets {A A') 
must be singlets. The inner product 'D-'D' of two doublets 
is also a singlet and invariant under any permutation 
because the representation matrices M G {M 12 , M 123 } 
are orthogonal: 

■ (MV) =V-V . (25) 

Therefore, there are three possibilities for constructing 
singlets in the product space: 

55', AA', V-V' = aa'Ass'. (26) 


5 Vj A 


Antisinglets are obtained from 



The multiplets do not mix under any permutation and 
thereby generate three invariant subspaces: 



( 21 ) 


® Our notation is based on Ref. [34]; see also [7, 35] for applications 
to the baryon Faddeev amplitude and the three-gluon vertex. 


SA!, AS', V hV-.= as'- sa', (27) 


where we defined an antisymmetric wedge product, and 
doublets are formed by 


SV, A{e'D'), _ as'A sa' 

S'V, A'{eV), * '“[aa'-ss'J’ 

where 


e = 





(28) 


(29) 


These relations are complete because they cover all 16 
possible combinations. One can verify them explicitly 
using Eq. (21): the resulting singlets stay invariant under 
permutations, the antisinglets pick up a minus sign for 
odd permutations, and the doublets transform under M 12 
and M 123 . 
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acb 



FIG. 3. Cayley graph for the permutation group S 3 , with 
permutations of positions rather than indices. 


B. Example: baryon flavor wave functions 


Let us illustrate these properties using the example of 
the well-known SU(3 )/ flavor wave functions for baryons. 
Their textbook construction can be somewhat involved 
but with the help of the relations above we can derive 
them in a few lines. 

The SU{3)f irreducible representations that are con¬ 
structed from the flavor vectors u, d, s in the fundamen¬ 
tal representation via 30303 = 108©8©1O differ 
by their symmetry, so we must classify them into simul¬ 
taneous irreducible representations of the permutation 
group S 3 . In this case we would like the group elements 
in Eq. (18) to interchange the positions of the objects 
a,b,c G {u, d, s} in a tensor product abc = a ® b ® c 
rather than the indices of / 123 . This requires slight mod- 
ihcations: in Fig. 2 we have to replace 


123 -G abc 213 —> bac 

231 —>■ bca 132 —> cba (30) 

312 —>■ cab 321 —> acb 

and hence P 13 <->■ P 23 in Eq. (18): 

1 , P12, 

P23Pl2=Pl23, Pl3=Pl2Pl23, (31) 

1^13 Pl2 = Pl23 ; P23 = P12 P123 ■ 



UUU uud ddu ddd 

UUS uds dds 

SSU ssd 

SSS 

5 

A++ A+ A° A- 

E+ SO S- 

so S- 

n- 

Pi 

P 2 

P ^ 

S+ SO E- 
AO 

SO S- 


A 


AO 




TABLE I. SU{3)f flavor wave functions for baryons. 


'spin-flavor ■ 


Ps ■ Pf 

(2,8) 

16 

^spin-flavor • 70 

Ss Sf 

(4,10) 

40 

Ps * Pf 

(2,8) 

16 

-^spin-flavor ■ 20 

Ss Pf 

Ps Sf 

(4,8) 
(2,10) 

32 

20 

Ps A Pf 

(2,8) 

16 

{ePs)Af 

(2,1) 

2 

Ss Af 

(4,1) 

4 





TABLE II. Spin-flavor wave functions from Eqs. (26-28). The 
brackets count the number of multiplets: the combination of 
two spin doublets Ps and eight flavor doublets Pf gives 16 
independent spin-flavor singlets, etc. 


together with P 2 =.4 = 0. Apart from overall normal¬ 
ization, S is the flavor wave function of the A+ and Vi 
that of the proton. Had we started from ddu instead of 
uud, we would have obtained the A*’ and the neutron (re¬ 
place u GG d above). The combination uuu returns only 
a singlet (A++), and from uds we get everything: S, A 
and two doublets. Taking into account all 10 combina¬ 
tions with different flavor content {uuu, ddd, sss, uud, 
uus, ddu, dds, ssu, ssd, uds), we arrive at the result in 
Table I: 

• ten symmetric singlets, which form the flavor decu- 
plet with A, E, S and H, 

• eight doublets that form the flavor octet, including 
proton, neutron, S, 2 and A, 

• and one antisymmetric singlet from uds, the flavor 
singlet for A. 


All other relations remain the same as long as the per¬ 
mutations are understood to act on positions rather than 
indices. 

To construct the flavor wave function for a baryon with 
flavor content uud, take a = b = u and c = d and apply 
Eq. (20) to the ‘permutation-group seed’ uud: 

= uud, ipp = 0 , = ±'0^ = - - - u. (32) 

Up to prefactors one obtains 


S = uud + udu + duu , 


Vi 


udu—duu 

— ^ (udu + duu — 2 uud) 


( 33 ) 


It is simple to extend the construction to SU (4) f, which 
produces 20 singlets, 20 doublets and four antisinglets 
(40404 = 205 © 20 m^ © 20 ms © 4^), or restrict it to 
SU{2)f (20 20 2 = 4s©2m^©2ms). 

The notation is especially convenient for constructing 
product wave functions. Take the full (Bethe-Salpeter) 
wave function of a baryon, 

pj = Dirac x Flavor x Color, (34) 

which must be totally antisymmetric under exchange of 
any two quarks. ’Dirac’ is here a shorthand for the full 
spatial-spin (or momentum-spin) contribution that trans¬ 
forms under the Poincare group. The color wave function 
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-^c = £ijk is totally antisymmetric. The flavor wave func¬ 
tions form multiplets 5/, 'Df and Af. Hence, the only 
combinations that produce a fully symmetric Dirac-flavor 
part are given by 


r {'DD-'Df)Ac (octet), 

Aotai = < {SDSf)Ac (decuplet), (35) 

[ {Ad A f) Ac (singlet). 

In the nonrelativistic quark model, the Dirac parts are 
the direct products of 0(3) orbital and SU{2) spin wave 
functions. The spin wave functions are identical to 
Eq. (33) and the first four columns in Table I if we re¬ 
place M by t and d by there are two doublets Vs with 
spin 5” = I and four singlets Sg with -S' = |. For orbital 
ground states {L — 0 ^ P = +, J = S) the orbital wave 
functions are spatially symmetric, so the only possible 
combinations are 

{Vs -Vf) Ac 

{SsSf)Ac 

whereas the flavor-singlet baryon does not appear in 
an orbital ground state. 

Combining the SU{2) spin multiplets Vg, Sg with the 
SU{3)f flavor multiplets Vf,Sf and Af yields the stan¬ 
dard SU{ 6 ) quark-model classification. The product 
wave functions can be read off directly from Table II: 
there are 56 singlets from Eq. (26) which are relevant for 
orbital ground states, 20 antisinglets from Eq. (27), and 
70 doublets from Eq. (28). 



(J = octet), 

(J = 1^, decuplet) 


IV. PERMUTATION GROUP Si 
A. Multiplets 

We will now generalize our notation from the last sec¬ 
tion to the permutation group S 4 . The group consists of 
4! = 24 elements. Each permutation of an object /1234 
can be reconstructed from two group elements, a transpo¬ 
sition P 12 and a 4-cycle T’i 234 - The former interchanges 
the indices 1 <->■ 2 and the latter is a cyclic permutation 
I—:^2, 2—>-3, 3—>-4, 4—>-1. For example, one has 

P23 = (^234)^ P12 7^1234 P12 , 

P34 = P12 (H1234)" P12 (P1234)" Pl2- 

The Cayley graph in Fig. 4 represents the group manifold 
as a geometric structure made of squares and hexagons. 
The squares contain the elements that are connected to 
each other by 4-cycles, for example 

f ^1234 f ^1234 f P 1234 f /qQ^ 

/1234 -^ /2341 -^ /3412 -^ 74123 ■ looj 

In total there are six squares with four elements each, 
which are connected to each other by transpositions P 12 . 
The two permutation chains in Eq. (37) then correspond 



FIG. 4 . Cayley graph for the group Si. The elements in the 
vertical center (red) are the t—channel permutations, the ones 
at the top and bottom (blue) are the u—channel permutations, 
and the remaining ones are the s—channel permutations. 


to paths along the Cayley diagram. Instead of P 12 and 
P 1234 it is also common to span the group by the transpo¬ 
sitions P 12 , P 23 and P 34 , which leads to a similar Cayley 
graph. 

One can rearrange the 24 permutations in multiplets 
that transform under the irreducible representations of 
S 4 . The respective Young diagrams are 

5 7 -+ V, rr A 



As for the case of 53 , 5 is a singlet and A an antisinglet, 
and the doublets Vj {j = 1,2) form a two-dimensional 
irreducible subspace. The triplets and antitriplets 
7”~ {i = 1,2,3) are new: they transform under inequiv¬ 
alent irreducible representations and thereby form two 
different three-dimensional subspaces. We denote their 
elements by 





(39) 


For the explicit construction of the multiplets it is help¬ 
ful to group the 24 permutations of an element /1234 into 
the following three subclasses: 


/1234 


/2314 


/3124 

/3412 


/4132 


/l 342 

/2143 

/l 423 

/4213 

_ /4321 _ 


. /3241 _ 


- 1 

CO 

A 
_ 1 

























together with the transpositions 

Pl 2 f^^\ Pl 2 f^^\ Pl 2 f^^^. (40) 

We label the column vectors with * = 1,2,3 and denote 
the column index by /c = 1... 4. Together with the trans¬ 
positions Pi 2 , each column defines a closed path in Fig. 4 
that contains eight elements. 

A similar classification is that in terms of t—, s— 
and M—channel permutations: the t—channel permuta¬ 
tions /(l), are the eight elements in the ver¬ 

tical center of the Cayley graph where {12} and {34} 
are grouped, the m— channel permutations , Pi 2 
are those at the top and the bottom with {31}{24}, and 
the s—channel permutations P 12 contain the re¬ 

maining eight elements with {23} {14}. Going from the 
t—channel to the s— and rt—channels amounts to 


/1234 


P 12 P 23 


> /2314, /1234 


P 23 Pl2 


> /3 


124 • 


(41) 


Now take the sum of the entries in each column and 
(anti-)symmetrize with P± = 1 ± P 12 : 


^t = p±J 2 f> 


(i) 


i = 1,2,3. 


(42) 


fc=i 


The resulting six combinations constitute the singlet, an¬ 
tisinglet and the two doublets which have the same form 
as the multiplets in S 3 : 

S = (V’l +'02 + 03)’’’ , 

-^ = (01 +02 + 03)~ , 


Vi 

V2 


(02 - 03) 

(02 +03 - 20l) + 

^ (02+03 - 20l)" 
(02 - 03)^ 


(43) 


The remaining 18 elements are reserved for the triplets 


and antitriplets. If 

we generalize Eq. (42) to 


= Px 

" 0 ^ 

1 

+ 

1 


hi = Px 

(/f ^ - /f - + /f) , 

(44) 

O 

II 

+ 

1 

1 


with A = ±, and define 


(0i)i = Xa^, 
(02)^ =Xbl 
(03)i = Ac^, 

(002 =-a2> (003 =a0 
( 02)2 = -c0 ( 02)3 = C 2 , 

( 03)2 = -hi, ( 03)3 = hi, 

(45) 


then we can cast the triplets and antitriplets in a common 
form: 




\Jl (01 + 02 + 03)f 
^ (02 + 03 - 20l)f 
(02 - 03)^ 


( 46 ) 


8 

The transformation laws for any permutation can be 
reconstructed from 



5 

-A 

M}2 Pj 


0 

1 

y 


S 

-A 


1^1234 

^{234 



(47) 


with the two-dimensional representation matrices 

Mi2=( I 5), ^1234 = ^^3 (48) 


and the three-dimensional matrices 


/ 1 0 0 

Hi2 = 0 1 0 

\0 0 -1 


^ / 2 4^2 0 \ 

H1234 = -X -2^2 1 3^3 

^ V 2^6 -V3 3 J 


(49) 


The representation matrices for the differ by a mi¬ 
nus sign from those of the 7}"'', and it is not possible to 
rearrange the entries {rt, u, ic} to obtain a common trans¬ 
formation law; i.e., the triplets and antitriplets transform 
under inequivalent representations.^ 


B. Product representations 

The next step is to find all possible products of ele¬ 
mentary multiplets {S, A, Vj and 7}^) that also satisfy 
the transformation laws in Eq. (47): 

(i) Singlets are obtained from the combinations 

55, AA, 

where (•) is the usual dot product for vectors: 

T) ■ V' := aa' + ss ', 

T -T' := uu' + vv' + ww'. 


(50) 

(51) 


^ An equivalent form of the representation matrices can be found 
in Table 3.2 of Ref. [361. They follow if one writes instead of 
Eq. (39): 



U 

+ 

W 

- 


-0 _ 

— V 

- r- = 

V 

, Vj = 

S 

a 


w 

i 

u 

i 



however at the price that Eqs. (46-47) and the notation for the 
product representations become less compact. The group charac¬ 
ters (the traces over the representation matrices) are necessarily 
the same. 
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The singlet property follows from the orthogonality of and there is a similar relation for the combination of 
the representation matrices. Similarly, the combinations triplet and doublet: 


5^, V.AVj, T^-T^ (52) 

produce antisinglets, where we used the antisymmetric 
product V AV = as' — sa'. 

(ii) Doublets can be constructed in various ways. S V 
and A {eV) are doublets; e was defined in Eq. (29). The 
combination of two doublets and two (anti-)triplets also 
produces doublets: 

T^*T^, e{T^*T^), (53) 


Ti • {T 2 AV)=T 2 - (Ti a V), 

Ti • {T 2 VV)=T 2 - {Ti V V). 

(64) 

Finally, we mention the identities 


{Vi A Vj) Vk = 0, 

(65) 

{ijk} 


E i-^f[T^-iT,ATk)]Ti = Q, 

(66) 

{ijkl} 



where the (*) operation has the form 


V*V' 

T*r 


as' + sa' 
aa' — ss' ’ 

vw' + wv' + '/2 [uw' + wu') 
ww' — vv' + {uv' + vu') 


(54) 

(55) 


(hi) Triplets T“'' are obtained from 


5r+, AT-, 


r+ vx>, 

T- AV, 


r± V , 

AT^ 


and antitriplets T from 


ST-, AT+, 


T+ AV, 

T- vv 


A , 

V . 


(56) 


(57) 


We defined the ‘wedge’ and ‘vee’ products for the triplets 
and antitriplets as 


TAT' 

TvT 


vw' — wv' 
wu' — uw' 
uv' — vu' 


vv' + ww' — 2uu' 
uv' + vu' + {vv' — ww') , 
uw' + wu' — {vw' + wv') 


(58) 

(59) 


and for the combination of (anti-)triplet and doublet: 


TAV 


T VV 


va — ws 

ua — ^ {va + ws) 
—us — ^ {vs — wa) 

vs + wa 

us — ^ {vs — wa) 
ua+ ^ {va + ws) 


(60) 


(61) 


The wedge product for the triplets is the usual vector 
product; as such it satisfies the identity 


Ti ■ {Tj A Tk) = Tj ■ {Tk A Ti), (62) 


where the sums are over cyclic permutations and P = ± 
for even/odd permutations. They are useful for expand¬ 
ing a doublet in a basis defined by two other doublets, or 
a triplet in a basis defined by three triplets. 

The list defined by Eqs. (50-57) is exhaustive. Check¬ 
ing the transformation properties of these quantities 
by hand can become tedious, but they are simple 
to implement in a computer algebra system such as 
Mathematica [37]. For example, suppose we want to find 
all possible products of one triplet T~^ with itself. The 
resulting six combinations u^, uv, uw, vw can 

be rearranged into a symmetric singlet, a doublet, and a 
triplet: 


T+ • T+ = -b -b , 


r+*r+ 


T+VT+ 


2 w {v + V2 u) 
w'^ — v'^ + 2\j2 uv 
-b — 2it^ 
2 uv + 'J2 (u^ — u>^) 
2 w {u — 'J2 v) 


(67) 


It is not possible to form an antisinglet, and the an¬ 
titriplet vanishes because T^ A T^ = 0. 

Table III collects all possible multiplets that can be 
systematically constructed from a doublet V and a triplet 
T^ =T up to cubic terms. This will become relevant in 
Sec. V because the six Lorentz invariants in a four-point 
function can be arranged into a singlet, a doublet and 
a triplet, which can generate further product multiplets 
with higher mass dimensions. 


C. Example: four-gluon vertex 

As an application, let us recast the color factors of the 
four-gluon vertex [38] in the multiplet notation. There 
are four types of color structures that can be used as seed 
elements for obtaining all further multiplets: 


where {ijk} is a cyclic permutation of {123}. It turns out 
that the identity also holds for the vee and star products. 


T-{TjVTk) = Tj-{TkVTi), 
Vi ■ {Vj * Vk) = Vj ■ {Vk * Vi), 


Aabcd — ^ab ^cd , 

Tabcd — fabe fade , 
^abcd — dabe ^cde , 
Vabed fabe ^ede ■ 


( 68 ) 
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(2,0) 

(1,1) 

(0,2) 

( 3 , 0 ) 

(2,1) 

(1,2) 

( 0 , 3 ) 

Singlet 

v-v 


T-T 

V-{P*P) 


V-{T*T) 

r- (TVT) 

Doublet 

V *7) 


T*T 

{V-V)V 


V*{T*T) 

{T-T)V 


Triplet 


T VT> 

rvr 


Tv {V*V) 
{V-V)T 

Tv {TVV) 

Tv (TvT) 
(T ■ T) r 

Antitriplet 


Tad 



Ta {V*V) 

TA (TVD) 

Ta (TvT) 

Antisinglet 




A (T> * T>) 


Va{T*T) 



TABLE III. Products of a doublet T> and a triplet T according to the rules ( 50 - 57 ) up to cubic terms. The brackets in the 
top row count the number of multiplets; e.g., ( 2 , 1 ) refers to two doublets and one triplet in the product. The elements in gray 
are trivially obtained by multiplying T> and T with singlets. 


Here, fate and dabc are the antisymmetric and symmetric 
structure constants of SU{N), respectively, 

/afec = Tr ([td, t^] tc), (69) 

date = 2Tr({ta,t;,}tc) , 

and ta are the SU(N) generators in the fundamental rep¬ 
resentation. 

We specialize to SU{3), where ta = Aq/ 2 and Aa are 
the Gell-Mann matrices. For the seed Aatcd = Sabred, 
the three column vectors in Eq. (42) simply become 


=SabScdy, 

/(") =Sbe6ady, V= } 

=5,aSbdy, Ll 


(70) 


due to the antisymmetry of the structure constants. The 
resulting antisinglet vanishes because of the Jacobi iden¬ 
tity, 

I A{B) = fabe fede + fbee fade + fcae fbde = 0, (73) 


and the remaining doublet is 


il?2(i3) = 


{fbee fade “t“ fcae fbde ‘^fabe fede^ 
fbee fade fcae fbde 


The seed Cabed = dabe dede produces another singlet and 
doublet X>i, but they linearly depend on the ones above: 


5(C) = i 5(H), 

I?i(C) = -|l?i(H) + Wp2(H), 


(74) 


and the transposition P 12 has the effect 


which is a consequence of identities that can be found in 
Ref. [38]. Finally, the seed Dabcd = fabedede produces 
one independent antitriplet and nothing else: 


Since all entries in V are identical, Eq. (44) can only 
produce zeros and therefore the triplets and antitriplets 
vanish. The ift of Eq. (42) become 




ifabe dede F fbee dade F fcae dbde^ 


— 8 dab ded ; 

tpf = 0 , 

f^2 ~ ~ ^ {dbe dad F dca dbd) ; 

^ {dbc dad dca dbd) ; 


(71) 


and hence only the singlet and the doublet Vi survives: 

g 5(H) — dab ded F dbc dad F dca dbd^ 

1 T) /' 4 ^ dad dca dbd 

® ^ {dbedad + dcadbd —‘2dabdcd) 


(72) 


The derivation for the seed Babcd = fabe fede is identical 
except that the transpositions P 12 produce minus signs 


with Tf~ (D) = —Tf~ (D). Therefore, the four-gluon ver¬ 
tex in 517(3) has eight independent color structures in 
total (we attach a subscript c for color): 


5c = 5(H), 


pW =Vi{A), 

vP =V2{B), 


Z-=Tf{D). (75) 


Ignoring the dabc symbols, they reduce to five. 

Now suppose we want to combine them also with 
the Lorentz structures. The four-gluon vertex is Bose- 
symmetric and therefore the products of Lorentz tensors, 
color factors, and momentum-dependent dressing func¬ 
tions must form symmetric singlets. If we restrict our¬ 
selves to the subset of momentum-independent Lorentz 
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FIG. 5 . Kinematics in the photon four-point function. 


structures, then only one seed element 5^'' con¬ 

tributes. It generates a singlet Sl and a doublet 
(with subscript L for Lorentz) which have the same form 
as in Eq. (72). Consequently, there are three possible 
Lorentz-color singlets in the product space: 

Vl-V^P, Vl-V^K (76) 

The last element is the tree-level structure of the four- 
gluon vertex because it can be written as 

= fate fcde - < 5 ^^ 5 "'^) 

+ fbcefadeiS^^S'''^-S>^-'Sn (77) 

+ fcaefMeiS>^''SP'^-S‘'P6n- 

If we include Lorentz tensors with higher mass dimen¬ 
sion, or if we allow the momentum-dependent dressing 
functions to form multiplets other than singlets, there 
will be many more possible combinations. A more de¬ 
tailed discussion of the Lorentz tensor basis will follow in 
Sec. VI. 


A. Definitions 

A generic vector four-point amplitude has the form 

N 

.. .P 4 ) = E /»(■■•) ■ •P 4 ), ( 78 ) 

i=l 

which is illustrated in Fig. 5 for the special case of four 
photons. It depends on four momenta pi, p 2 , Pa, P 4 , of 
which only three are independent due to momentum con¬ 
servation Pi +P 2 +P 3 +P 4 = 0. Since the Lorentz indices 
are not relevant for the discussion in this section, the 
analysis applies not only to the hadronic LbL amplitude 
but also to the four-gluon vertex or any other symmet¬ 
ric four-point function. M has a decomposition into N 
tensor basis elements with N Lorentz-invariant dressing 
functions fi. They depend on 4 • 3/2 = 6 independent 
Lorentz invariants which we will specify below. 

Bose symmetry entails that all 24 permutations of A4 
(interchanges of momentum labels and Lorentz indices) 
are identical: 

7W'"'"'’‘"(P1,P2,P3,P4) = M''^P‘^{p2,Pl,P3,Pi) = ■■■ 

We can group them into three classes with eight permu¬ 
tations each, as discussed earlier: t-channel {12}{34}, 
s-channel {23}{14} and u-channel {31}{24}. Therefore, 
the simplest choice of independent momenta are the t—, 
s—, and u—channel Mandelstam momenta: 

k=Pi+P2 = -P3 - Pi 1 
p = P2 + P 3 = -Pi - P 4 , ( 79 ) 

q=P3+Pi = -P 2 - Pi , 

with the inverse relations 

Pi =\{q-P + kl), P 2 = -\{q- p- k), 

P3 = \{q + p-kl), P4 =-\{q+p + k). 

Going from the t—channel to the s— and u—channels 
then amounts to cyclic permutations 


V. PHASE SPACE IN FOUR-POINT 

t — 

channel : 

Pi, 

P 2 , 

P3, 

Pi 


P, 

9 , 

fc. 


FUNCTIONS 

s — 

channel : 

P 2 , 

P3, 

Pi, 

Pi 


9 , 

fc. 

P, 

( 81 ) 


u — 

channel : 

P3, 

Pi, 

P 2 , 

Pi 


fc. 

P, 

9 - 



We proceed by applying the permutation-group tech¬ 
nique to the phase space in four-point functions. As we 
discussed in Sec. II, in a Bose-symmetric and minimal ba¬ 
sis all form factors are free of kinematic singularities or 
zeros and their momentum evolution is governed purely 
by the dynamics of the system. A generic offshell four- 
point function has a rich phase space because its form 
factors depend on six Lorentz invariants. Therefore, the 
goal of this section is to isolate the relevant kinematic re¬ 
gions where dynamical singularities occur, and to identify 
possible ‘scaling variables’ that simplify the description 
because they carry the main momentum dependence. 


B. Phase space 

In the next step we would like to develop a geometri¬ 
cal understanding of the phase space, in particular of the 
spacelike (‘Euclidean’) region. In the context of the LbL 
scattering amplitude, this is the integration domain that 
contributes to the muon g — 2 value, whereas in the four- 
gluon vertex it is the domain that is mapped into itself in 
a Dyson-Schwinger equation. The fi in Eq. (78) depend 
on six Lorentz invariants that constitute the phase space: 


12 


from the momenta above one can form the three Man¬ 
delstam variables , q'^, and the angular variables 

oji = q ■ k, 0 J 2 = p ■ k, (jj 3 = p ■ q . (82) 

The Mandelstam variables are convenient for discussing 
two-photon intermediate states that appear in the LbL 
amplitude, as we shall see below. 

To expose the Bose symmetry of the phase space, we 
would like to arrange these variables into permutation- 
group multiplets. The simplest strategy is to first ar¬ 
range the four-momenta p, q, k themselves into multi¬ 
plets. If we use the t—channel momentum k = pi + p 2 
as permutation-group seed /1234 and follow the steps in 
Eqs. (42-46), the only nonvanishing multiplet turns out 
to be the triplet 


r+ = 


1 

2 


:^{p + q + k) 

^{p + q-2k) 


(83) 


These are just the three independent Jacobi momenta of 
the system. Applying the three product operations in 
Eq. (67) to T'*’, where each operation is simultaneously 
a scalar product of four-momenta, yields six Lorentz in¬ 
variants that form a singlet, a doublet and a triplet. The 
singlet is given by 


whereas the doublet is 


a 

1 

■ V3(q2-p2) ■ 

s 

“ 4 

_ p^ -1- q^ — 2k‘^ _ 


and the triplet has the form 



u 

1 

A 

-2 {oji - 1 - 012 + Wa) 

II 

V 

—-72 (wi -|- 012 — 2 + 13 ) 


w 

4 

Vq {oji - CLI 2 ) 


(84) 


(85) 


( 86 ) 


We pulled out factors of Sq to remove the mass dimension 
of the doublet and triplet variables, so that only 5o G 11+ 
carries a dimension. The advantage is that the doublet 
and triplet phase spaces {o, s} and {m, u, ui} can be dis¬ 
cussed independently of 5o (otherwise they would scale 
with 5o). 

In Sec. VI we will express the photon four-point func¬ 
tion in a fully symmetric basis, so that the corresponding 
dressing functions must be also symmetric and can only 
depend on singlet Lorentz invariants. The six lowest- 
dimensional singlets constructed from Sg, Vg and To are 
those in the first row of Table III: 


'Tl — Vg ■ Vg , 

S3 = To • To , 


S2 = TJq ■ i'P’o * T^o) , 
54 = Vg ■ {Tg * To) , 
Sg = To ■ (To V To), 


( 87 ) 


together with 5o itself. Since 5o has the lowest mass di¬ 
mension, it is reasonable to assume that it also carries the 
main momentum dependence, whereas the dependencies 
in the remaining variables are weaker. 

What is the domain for {a, s} and {m, v, to} in the 
spacelike region? To find out, we express the three Man¬ 
delstam momenta in hyperspherical coordinates (using 
Euclidean conventions): 



roi 


0 

II 

0 

0 

73 

II 

to 1 

0 

Vl - 


L 1 J 


Z 


T = 


- 0 

v ^l - \/l - 

\/l - y 


( 88 ) 


The radial variables and q^ are all real and posi¬ 

tive, whereas the variables z, z' and y are the cosines of 
hyperspherical angles: 

p^,q^,k'^ e E+, z,z',y G [-1,1]. (89) 

It is straightforward to derive their relations with the 
Lorentz-invariants in Eq. (82). It turns out that {a, s} 
form the interior of a triangle and {u, u,w} that of a 
tetrahedron, which we will discuss next. 


C. Doublet 


The doublet phase space is the Mandelstam plane 
spanned by the variables a and s defined in Eq. (85). 
It encodes the relations between p ^, q^ and k^ which are 
illustrated in Fig. 6. The spacelike region defined by 
Eq. (89) forms an equilateral triangle of side length 2-\/3, 
enclosed by the lines = 0, = 0 and = 0. The 

three corners 


a 


0 


[73] 


[-73] 

s 


-2 

5 

1 


1 


correspond to the limits q^ = p^ = 0, p^ = = 0 and 

q^ = = 0, respectively. Since the three momenta are 

independent and the variables p^, q^ and can take any 
value G E+, all points within the triangle contribute to 
the spacelike region.® The center of the triangle is the 
symmetric point where p^ = q^ = fc^. 

In practice it is convenient to parametrize the doublet 
by a radius and an angle: 


2To — ^0 



sinp 
— cos Lp 


(91) 


^ For a three-point function one obtains a similar doublet, namely 
a triangle in the variables and pg. However, since in 

that case only two momenta are independent, the phase space is 
restricted to the interior of a unit circle embedded in the triangle 
(indicated in Fig. 6), cf. Ref. [7]. 
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FIG. 6. Doublet phase space for the photon four-point func¬ 
tion in the (a, s) plane at a slice of fixed 5 . 


with (p G (0, 27r), chosen such that (p = 0 corresponds to 
the lower corner of the triangle. Note that the range or r 
is a function of the angle ip as spelled out below, Eq. (94). 
When expressed in r and <p, the Mandelstam variables are 
cyclically related by a rotation ip ^ p± 27r/3: 

p" = i|?{l + rc«{^ + 2i)), 

^(l+>’Cos{v=-^)), (92) 

2 4iSo /, N 

k = —— (1 -I-r cost/?). 

O 

The sum of the cosines vanishes, so that -\- k^ = 

4iSo from Eq. (84) is satisfied. The permutation-group 
invariants in Eq. (87) become 

= Sq , S 2 = Sq cos 3i^ . (93) 

The transformation of the doublet under any permu¬ 
tation of momenta can be reconstructed from the Cayley 
diagram in Fig. 4 and the representation matrices M 12 , 
M 1234 in Eq. (48). Permutations act on the angle ip only: 
P 12 induces a reflection ip —>■ —ip and P 1234 a rotation 
and reflection ip —>■ —ip — 27r/3. Therefore, permutations 
within the t, s or u channels reflect the triangle along 
any of the three axes (dashed lines) in Fig. 6 ; permuta¬ 
tions between the channels rotate it by an angle ±27r/3. 
The maximum value of the radius r (the boundary of the 
triangle) defines the boundary of the spacelike region: 

^max(¥’) = r ' . I 3 m iNt" ‘ (94) 

This is again a permutation-group invariant function 
with 1 < rniax(<p) < 2. The variable f = r/r^^-y^{ip) 
takes values in the interval [ 0 , 1 ] only. 


The Mandelstam triangle is convenient for discussing 
two-photon intermediate states that appear in the LbL 
amplitude. While the space-like region is free of singu¬ 
larities, it is still influenced by the time-like singularity 
structure. For example, it is well known that the pion 
poles in the two-photon channel produce a large (and pre¬ 
sumably even the largest) QCD contribution to the LbL 
amplitude [I]. Bose symmetry entails that two photons 
can only produce meson quantum numbers with even 
charge-conjugation parity: = 0 ^ (pseudoscalar), 

0++ (scalar), I++ (axialvector), I ^ (exotic vector), 
2++ (tensor), etc. Such poles will occur at p^ = —w?, 
= —m? or = —m?, where m is the respective pole 
mass (which is complex if the state carries a width). In 
terms of the radius this translates to 

f=I+—>1. (95) 

40o 

Hence, each real pole constitutes another triangle that 
encompasses the one in Fig. 6. They should have an 
impact on the spacelike region and induce a rise in the 
dressing functions towards the boundary f = 1. Because 
So appears in the denominator of Eq. (95), the sensi¬ 
tivity to timelike physics will be generally stronger at 
larger values of So, i.e., in the UV region, because the 
poles approach the boundary of the triangle, whereas for 
iSo —>■ 0 they move infinitely far away. A similar behavior 
is known from the three-gluon vertex, where the angular 
dependence is sizeable in the UV but negligible in the 
infrared [7]. 


D. Triplet 


Whereas the Mandelstam triangle encodes the two- 
photon intermediate states, the triplet variables {u, v, re} 
are related to the photon virtualities Xi = p], i = 1 ... A, 
and thereby encode the vector-meson poles. To see this, 
one derives the following relations between the Xi, So and 
Wi from Eq. (80): 


Xi = 5o + 4 (wi -UJ2- W 3 ), 
X2=So+l (-Wi +UJ2- UJ3) , 
0:3 = 5 o + 4 (-Wi - W2 + UJ3), 

X4 = Sq + ^ (wi + UI2 + UJ3) . 

In terms of the Xi the singlet takes the form 


5o = 


Xi + X2 + X3 + X4 


and the triplet is given by 


To = - 


Xi + X2 + X3 — 3X4 

— V2 {xi + X2 — 2X3) 
VQ{xi - X2) 


(96) 


(97) 


(98) 


If we allow the Xi to take any spacelike value S 1R,+ , 
the resulting phase space forms a tetrahedron which is 
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FIG. 7. Triplet phase space for the photon four-point func¬ 
tion in the (u, v, w) space. 


shown in Fig. 7. Its four corners To = SgCi, i = 1... 4, 
are defined by the kinematic limits where Xi = 4iSo and 
the other three virtualities Xj^i vanish: 



1 


1 


1 


■ -3' 

C^ = 

-V2 

1 



2^/2 

1 

0 


V6 


-v^ 


0 


0 


(99) 

The soft kinematic limits where only one of the photons is 
real {xi = 0 ) dehne the four faces of the tetrahedron; i.e., 
they are the planes spanned by the points Cj^i. These 
are the points which are relevant for {g — 2 )^, since the 
external source of the magnetic held is represented by 
an on-shell photon [1]. The center points of the four 
faces Ai = —CilZ lie on the opposite sides of the corners 
Ci- The center of the tetrahedron is the symmetric point 
where all Wi vanish and consequently all Xi = Sq are the 
same. 

Similarly to Eq. (91), we parametrize the triplet in 
spherical coordinates, 


For the Wi this choice entails 

wi = R (^cos9 + V 2 sin 0 cos [4> + , 

W 2 =-^ i? ^cos0-I-sin0 cos ((/) — ^) j , (101) 

W 3 = R ^cos 9 + \f2 sin 9 cos (/>^ , 

which are again cyclically related by a rotation ^ ^ ± 

2tt jZ. The remaining invariants from Eq. (87) then take 
the form 


S^ = SlR\ ( 102 ) 

^4 = Sq rR? [sin^ 9 cos(2(/) -|- (^) -|- V2 sin 20 cos(^ — ip )\, 
1 S 5 = Sq R^ cos 30 -I- cos^ 0 -I- V 2 sin^ 0 cos 3(j )). 

The representation matrix H 12 for the permutation P 12 
induces only a rehection (j) —>■ —(j), whereas H 1234 is more 
complicated to express in these variables because it ex¬ 
changes the faces of the tetrahedron. Still, going to the 
s-channel (P 12 T 23 ) or rt-channel (P 23 P12) is simple be¬ 
cause it only induces a rotation (j) —>■ ^ ± 27r/3. This 
is quite helpful in practical calculations of the four-point 
function when the dressing functions are expressed in the 
variables So, r, ip, R, 0, (f>. Once the sum of all t—channel 
diagrams has been obtained, the remaining permutations 
in Eq. (81) amount to nothing more than a simultaneous 
rotation of the angles p and (j) by ±27r/3 each. 

A complication comes from the fact that not the whole 
interior of the tetrahedron contributes to the spacelike 
region. Only three of the momenta pi are independent, 
which leads to a restriction on the phase space.® For 
given doublet variables r and p, the actual triplet domain 
is a complicated geometric object contained within the 
tetrahedron. It is visualized in Fig. 8 for five different 
kinematic conhgurations. Its surface is defined by the 
permutation-group invariant condition^ 

p^q^k"^ = ujIp^ -\- 0 J 2 + oj'lk'^ — 2 wi UJ 2 0 J 3 , (103) 


which leads to a cubic equation for the radius R. Its 
solution is: 


Rmayi(r,p,9,(j)) = 


1 -I- 2sgn(X) cos 


$ 


$ = arccos |X| -|- 27r [0(A) 0(—A) — 0(A)]. 


(104) 


Here, 0 is the unit step function and the remaining aux¬ 
iliary quantities are the invariants 


A = ^^ 
253 A 3 ’ 


B = l + 


Si 

2S^R^ 


(105) 



U 


— cos 9 

% — So 

V 

— So R 

sin 0 cos (j) 


w 


— sin 0 sin (jj 


0 e ( 0 , 7 r) , 

(/) G (0,27r) . 

( 100 ) 


® This is similar to the restriction for the doublet phase space in a 
three-point function, see footnote 5. 

^ Expressed by the variables in Eq. (88) this corresponds to either 
= 1 or = 1. 
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r = 0.9 r = 0.5 

(p = n (p = n 




r = Q 





r = 1 r = 1.6 

(p = 0 (p = 0 


FIG. 8. Different kinematic points in the doublet triangle (bottom row) and the corresponding spacelike region within the 
triplet tetrahedron (top row). The alignments are the same as in Figs. 6 and 7. A rotation ip ^ ip ± ^ will also rotate the 
tetrahedron. The figure in the center (r = 0) contains both the central limit {R = 0) and the limit of three equal photon 
momenta {R = 2, the corners of the spacelike volume). 


with domains Ad [—1,1], B G [0, 3] and 


X = 1- 


2A^D 


D = ( 106 ) 


with domain D G [0,1]. Consequently, i?max is also in¬ 
variant under permutations. It is restricted to the in¬ 
terval 0 < i?niax ^ 2 and never reaches the corners of 
the tetrahedron with radius i? = 3. In other words, 
three photon virtualities can never vanish simultaneously 
in the domain that is integrated over in the g — 2 cal¬ 
culation. The variable R = R/R„ia.yi(j',ip,9,4>) is again 
permutation-group invariant and its domain is the inter¬ 
val [0,1]. 

Since the triplet phase space encodes the relations be¬ 
tween the photon virtualities, additional momentum de¬ 
pendencies will come from vector-meson poles at time¬ 
like virtualities Xi = —rn^. Microscopically, each photon 
in the hadronic part of the LbL amplitude couples to a 
quark via a dressed quark-photon vertex, which automat¬ 
ically contains all vector-meson poles. The sequence of 
vector-meson poles corresponds to a sequence of tetra- 
hedra that encompass the one in Fig. 7. In analogy 
to the doublet, the dressing functions in the spacelike 
domain should therefore rise in magnitude towards the 
spacelike boundary. Fig. 8 shows that the sensitivity to 
vector-meson poles should be weakest in the corners of 
the doublet triangle (where the triplet volume shrinks 
to the point R = 0, rightmost figure) and strongest at 
the centers of its sides (where the triplet extends to the 
edges of the tetrahedron at i? = 2, leftmost figure). Gen¬ 
erally, while the spacelike interior of the triangle and the 


tetrahedron is free of singularities, one- and two-photon 
singularities will therefore influence the behavior of the 
dressing functions from the timelike domain. 


E. Special momentum configurations 

The main practical goal of the permutation-group anal¬ 
ysis was to facilitate the discussion of the phase space. 
Once a symmetric tensor basis is employed, the dressing 
functions can only depend on the symmetric combina¬ 
tions in Eq. (87). If the basis is minimal, it fully absorbs 
the kinematic part of the momentum dependence. Hence, 
the momentum evolution of the form factors should be 
mainly governed by the scaling variable Sq with lowest 
mass dimension: fi{So ■. - S^) « fi{So), and therefore be¬ 
come simple. However, angular dependencies induced by 
timelike singularities will become relevant towards the 
boundaries of the triangle and tetrahedron. 

We conclude this section by discussing some special 
momentum configurations, some of which have been also 
studied in the context of the four-gluon vertex [12], Ex¬ 
pressed in our permutation-group variables, they corre¬ 
spond to the following limits: 

(i) Uniform soft limit {Sq = 0): all momenta vanish 
simultaneously, p = q = k = Q, and consequently also all 
Lorentz invariants are zero. 

(ii) The soft-photon limit is the relevant limit for 
the g — 2 problem: one external momentum vanishes, e.g. 
P4 = 0 and therefore p + q + k = Q. The variables q^ 
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FIG. 9. Phase space that is relevant for g — 2. 

and are still independent but 

wi = 5 - k^), 

u;, = -l(p^-q^ + k^), (107) 

In terms of doublet and triplet variables, Sq, a and s 
remain independent whereas 

u = I, V = —V^s, w = —\f2a. 

This is the lower face of the tetrahedron, whose remaining 
variables v and w are now proportional to those in the 
Mandelstam triangle s and a. The triplet variables from 
Eq. (100) become 

i ?2 = l + 2 r^ cos 6 » =- ^ o = 

Vl + 2r^ 

Since i? = i?max and therefore i? = 1, the combination 
with the upper limit for i?max from Eq. (104) restricts 
the Mandelstam triangle to the unit circle: r < 1. The 
domain is visualized in Fig. 9. 

(hi) Central limit: the Mandelstam momenta have 
the same length and are orthogonal to each other: 

= I 5o, a;i=W 2 =W 3 = 0 . (108) 

Therefore only 5o 7 ^ 0, whereas a = s = u = v = w = 0. 
This is the center of the triangle and the tetrahedron. 
The doublet and triplet radii vanish, r = i? = 0, and all 
photon virtualities are equal: Xi = Sq. 

(iv) Three equal momenta: pi = P 2 = Pa = —pi/'i: 
hence all Mandelstam momenta are identical {p = q = k) 
as well as the Lorentz-invariants: 

p^ = g^ = = wi = a ;2 = wa = I 5o . (109) 

It entails a = s = v = w = 0 but u = —2; the photon 
virtualities are Xi = X 2 = X 3 = ^ Sq and X 4 = SiSq. 

(v) Two equal momentum pairs: for example. 
Pi = P 2 = —P 3 = —P 4 - As a consequence, p = g = 0 and 
therefore a = u = v = w = 0 and s = — 2 , which is the 
limiting case for the right panel in Fig. 8 . 


(vi) Two opposite momenta: e.g. pi = —p 2 and 
therefore pa = —P 4 , so one of the Mandelstam momenta 
vanishes (fc = 0). There are now three independent vari¬ 
ables: p^, g^ and wa = p • g, with k'^ = loi = 1 x 2 = 0. 
Equivalently, s = 1, la = 0 and v = —y/2u, which is the 
limiting case for the left panel in Fig. 8. 

(vii) Two equal momenta: e.g. pi = P 2 and there¬ 

fore also p = g. Also here three independent variables 
remain: p^ = = wa, and p ■ k = loi = UJ 2 , which 

corresponds toa = w = 0, '!; = (n-|-s-|- 2 )l\f 2 . 


VI. TENSOR BASIS 

Having discussed the phase space, we will now proceed 
with the tensor basis construction for the photon four- 
point function. It should respect the symmetry prop¬ 
erties of the permutation group 54 and the transversal- 
ity and analyticity constraints. Gauge invariance entails 
that the LbL amplitude is transverse, and analyticity im¬ 
plies that it must be at least quartic in the photon mo¬ 
menta. The form factors in such a basis should be free of 
kinematic singularities but also free of kinematic depen¬ 
dencies and kinematic zeros. If the tensor basis is made 
of permutation-group singlets, the form factors must be 
symmetric and can only depend on the singlet variables 
in Eq. (87), so that their momentum dependence should 
become simple. 

In Sec. II we worked out the simpler case of a scalar 
two-photon current which will serve as a template in what 
follows. However, we will encounter additional subtleties 
that do not appear in the two-photon example. In this 
respect, a closer analogue is nucleon Compton scatter¬ 
ing, whose tensor basis was worked out by Bardeen and 
Tung [23] and Tarrach [24]. In that case the procedure 
is as follows: (i) Write down all possible 34 tensor struc¬ 
tures that are compatible with Lorentz covariance. It 
turns out that two of them are redundant and can be 
eliminated, thus leaving a 32-dimensional type-I basis, 
(ii) Apply gauge invariance to arrive at a transverse ba¬ 
sis with 18 elements. Tarrach noted that the resulting 
basis (Eq. (12) in Ref. [24]) was not minimal, i.e., the 
basis elements were still kinematically dependent, and 
suggested to replace three tensors by alternative ones 
in a certain kinematic limit (Eq. (13) therein). It was 
later realized by the Mainz group [26] that this would 
not have been necessary had one started directly with 
the 18-dimensional basis made from Tarrach’s three al¬ 
ternative tensors, because that basis is minimal. The 
reason behind this becomes only clear when the basis el¬ 
ements are cast into permutation-group singlets, because 
the corresponding three new tensors have lower powers 
in the photon momenta than the original ones [39].® 


What complicates the power counting in Compton scattering is 
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We will make similar observations also for the photon 
four-point function. In Sec. VIA we construct a generic 
type-I basis with 138 elements, but it will turn out that 
two of them are redundant and can be eliminated. In 
Sec. VIB, after applying the constraints from gauge in¬ 
variance, we will arrive at 41 transverse tensor elements. 

The permutation-group technique will be an essential 
ingredient in determining the minimality of such bases. 
The criterion is that the singlet basis elements should de¬ 
pend on the lowest possible powers of photon momenta, 
which translates into the lowest mass dimensions. In the 
Compton scattering example, invariance under charge 
conjugation and photon crossing leads to an S 2 x S 2 sym¬ 
metry, and constructing a singlet basis is simple because 
it is sufficient to multiply the (anti-)symmetrized basis 
elements with appropriate momentum prefactors. For 
the photon-four point function the situation is more dif¬ 
ficult because of the more complicated structure of the 
group S 4 : one has to cast the tensor basis elements into 
irreducible representations of S 4 , combine them into sin¬ 
glets, and ensure that the resulting basis elements have 
the lowest possible mass dimension. 

In any case, we proceed under the assumption that 
it should be possible to find a 41-dimensional minimal 
basis without the need for swapping basis elements in 
different kinematic limits. Since the problem is quite 
formidable we have not yet succeeded in this goal, but 
we will describe the necessary steps in the following. 

A. Type-I basis 

To begin with, we prove that the LbL amplitude (and 
also the four-gluon vertex) has indeed 136 independent 
Lorentz tensors and not 138. The system depends on 
three independent momenta qi (i = 1,2,3); for exam¬ 
ple, we can work with the three Mandelstam momenta 
in Eq. (79). A naive counting of all possible combina¬ 
tions of Kronecker deltas and four-momenta yields 138 
elements: 

• ^ 3 permutations ^ 3 elements, 

• Qj 3^ X 6 permutations = 54 elements, 

• Qi Qj Qk Qi ^ 3"^ = 81 elements. 

The list is, however, redundant: for four- and higher 
n-point amplitudes the dimensionality of spacetime re¬ 
stricts the number of independent basis elements. This 
can be understood from a simple argument. Suppose we 
orthogonalize the three independent momenta to obtain 
three unit vectors rii, i = 1,2,3 that are transverse to 


that the photon momentum powers are no longer equal to the 
mass dimension: the amplitude also depends on the average nu¬ 
cleon momentum, and the contraction with nucleon spinors can 
additionally raise these powers due to Gordon identities. 


each other. From three vectors one can construct an ax- 
ialvector niU^nl with opposite parity. A 

complete, orthonormal tensor basis follows from collect¬ 
ing all possible combinations of 

• v" ^1 element, 

• v'^ Uj => 3^ X 6 = 54 elements (HO) 

• n'' n^n'^ ^ 3^ = 81 elements, 

where v must appear in pairs to ensure the correct parity. 
This yields 136 tensor structures instead of 138. The 
list is already complete since the Kronecker delta can be 
written as the linear combination 

3 

v’'+ n''. (Ill) 

i=l 

In a frame where v, ni, 712 and are the unit vectors in 
IR 4 this is obviously true; since the equation is Lorentz- 
covariant it holds in any frame. Therefore, the element 
permutations are already contained in the 
basis ( 110 ) and do not generate additional structures. 

We emphasize that this reduction is a consequence of 
dimensionality and has nothing to do with gauge invari¬ 
ance (which we will exploit in Sec. VIB) or Bose symme¬ 
try (which does not reduce the number of independent 
tensors). One can repeat the exercise for higher n—point 
functions: also in that case at most three of the (n — 1 ) 
independent momenta can appear in their bases, which 
greatly reduces their dimensions. 

The same argument also gives the correct number of 
transverse elements. When applying transverse projec¬ 
tors, no elements with the same label and index survive: 
Pi is longitudinal, and so are P 2 , p^ and pf. At least 
for counting purposes, removing these four objects from 
the basis is formally equivalent to crossing off any of the 
vectors Ui from the list above, so that * = 1,2 only. This 
leads to 1 -I- 2^ X 6 -I- 2^ = 41 transverse elements. 

We note that the basis in Eq. (110) is already orthonor¬ 
mal because all the vectors and axialvectors it contains 
are normalized and orthogonal to each other. Hence, it 
defines the simplest possible complete tensor basis for 
a vector four-point function. On the other hand, nei¬ 
ther Bose symmetry nor gauge invariance or analyticity 
are implemented at this point. In the aforementioned 
Lorentz frame all basis elements are non-singular because 
they consist of pure numbers, but the corresponding 
(Lorentz-invariant) dressing functions will exhibit kine¬ 
matic dependencies and zeros. 

Therefore, we should first construct the analogue of 
Eq. (4): a 136-dimensional type-I tensor basis that is 
free of kinematic singularities and made of permutation- 
group singlets. It will serve as our starting point for 
working out the transversality constraints. The simplest 
construction principle is to use the Mandelstam momenta 
p, q and k together with the permutation technique from 
Sec. IV A, i.e., to write down the maximum number of 



18 


seed elements and work out their permutations. The out¬ 
come of this procedure is collected in Table IV. For ex¬ 
ample, the permutations of the seed element 5^'^SP'^ are 

SPPS'''^. ( 112 ) 

With Eqs. (42-46) they can be arranged into a singlet S 
and a doublet Vi whose structure is analogous to those 
in Eq. (72); all other multiplets vanish. Likewise, the 
seed SP'' kPk'^ generates the six permutations 

SP''kPk'^, 5'^PpPp'^, SPPq'^q'^, 

dP'^kPk’', SP'^p^pP, 5'"^qPqP ^ ^ 

which produce a singlet, a doublet and a triplet T{^. 
The remaining Lorentz tensors with one Kronecker delta 
and two identical momenta come from the seed S^'^pPp'^ 
(or equivalently S^'^qPq^) which has 12 permutations: 

SP'^pPp'^, 5''PqPq’^, 5PPk'^k^, 

SP'^pPp'', SP'^q^'qP, 5'"^kPkP 

(114) 

6P''qPq'^, S'^PkPk'^, 5PPp''p’^, 

SP’^qPq’', 6P’^k''kP, 5''‘^pPpP. 

One proceeds along these lines until the list is complete. 
The resulting 136 Lorentz tensors in Table IV are ar¬ 
ranged with increasing mass dimension n\ there are 3, 
54 and 79 elements for n = 0, 2,4, respectively. Because 
each seed produces only one symmetric singlet and there 
are 11 singlets in total, this is also the minimum number 
of independent tensor elements in the four-point func¬ 
tion: all 136 tensors can be reconstructed from those 
eleven through permutations. 

In principle, for n = 4 there would be 81 independent 
elements but here the spacetime restriction discussed 
above comes into effect: the basis saturates with 136 
elements and adding two more produces linear depen¬ 
dencies. Which ones to remove is not arbitrary because 
unfortunate choices can produce kinematic singularities 
in the dressing functions already at this stage. If we label 
the eleven seed elements in Table (IV) by 1... 11 from 
top to bottom, then one can show that the basis element 
Vi (8) — which is the one in brackets in the table — is a 
linear combination of the multiplets 

5(2)-5(3), 5(7)-5(8), 

I?i(l), I?i(2), I?i(3), I?2(3), I?i(7), I?2(8), 

r+(4) - 7;+(5), r+(9) - r 3 +(ii), ( 115 ) 

r2-(5), r 3 -(ii), 

y4(3), Ai8). 

The coefficients are rather lengthy but they respect the 
doublet construction rules in Eq. (53). It turns out that 
all coefficients share the denominator ^ (r^ — 16), where 
r is the doublet radius defined in Eq. (91). Hence, unless 
r = 4 (which is never reached in practice because r < 2 


n 

Seed 

# 

Multiple! type 

0 

SP'^SP'^ 

3 

5, Vi 

2 

&P'' kP k” 

6 

5, I?i, T+ 


SP'^pPp'^ 

12 

5, Vi, ©2, T^, A 


SP^'pP g'" 

12 

5, Vi, 7)+, 7^* 


SP''pP k'^ 

24 

5, ©1, ©2, r±, r2±, Tt, A 

4 

ppp''pPp‘’ 

3 

5, Vi 


pp p" qP 

6 

iS, 7^1, 1\ 


pp p" kP k’’ 

10 

S, (Di,) ©2, r±, .4 


pP q-' kP k” 

12 

5, T>1, 7)+, 


pPp'^ pP fe" 

24 

5, ©1, ©2, T±, r2±, Tt: .4 


pp p'' qP k’’ 

24 

5, T>i, T>2, T^, T^, T^, a 


TABLE IV. 136-dimensional tensor basis for the vector 
four-point function, where gauge invariance is not yet imple¬ 
mented. The doublet in brackets is linearly dependent due to 
the spacetime restriction discussed in the text; its inclusion 
would lead to 138 instead of 136 tensor structures. 

in the spacelike domain) the element T>i(8) is linearly 
dependent and can be removed.® 

In principle we still need to recast the tensor struc¬ 
tures in Table IV into permutation-group singlets. Since 
the procedure is the same for the type-I basis and the 
transverse basis that we will derive next, we integrate 
the discussion into the following subsection. 

B. Transverse tensor basis 

The remaining task is to work out the consequences 
of electromagnetic gauge invariance. We start from the 
expression (78) for the LbL amplitude, 

136 

MP'^P^p, q,k) = Yl 9 ’ > ( 116 ) 

where the {p^q^k) are the basis elements from Ta¬ 
ble (IV) or, alternatively, 136 singlets constructed from 
them. The transversality conditions have the form 

... = ^ (117) 


^ This is not entirely satisfactory but sufficient for our present pur¬ 
poses. Ideally it should be possible to remove two elements in 
arbitrary kinematics, as it can be done for the Compton scatter¬ 
ing amplitude [24]. 
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and reduce the basis to a subset of 41 transverse tensors. 
Transversality and analyticity require these tensors to be 
proportional to at least four powers in the photon mo¬ 
menta. Instead of Eq. (117), one can equivalently work 
out the Bose-symmetric condition 

Tf“ rp Tf , (118) 

where the Tp = p‘'/pi are the transverse projec¬ 

tors with respect to each photon momentum. This leads 
to relations between the dressing functions; if we denote 
the independent functions by /i and the dependent ones 
by they take the form 

9i = 5i(/i! ■ ■ ■ / 41 ), 

: (119) 


decomposition ( 120 ) provides a convenient filter that al¬ 
lows one to quantify such gauge violations and, if possi¬ 
ble, remove them to arrive at physically meaningful pre¬ 
dictions. 

While the procedure outlined here is at least in prin¬ 
ciple straightforward, it is almost impossible to perform 
by hand because of the sheer length of the expressions 
involved. Hence, we take the alternative route that we 
advertised in Sec. II, which is the essence of Tarrach’s 
procedure [24]: construct tensors with lowest possible 
mass dimensions that are automatically free of kinematic 
singularities. The mass dimension must be even because 
the four-point function has positive parity. By working 
out the permutations of these tensors we can construct 
a linearly independent, complete basis made of 41 ele¬ 
ments. 

To this end we employ 


995 — 595(/i, ■ ■ • /41) ■ 

They must be solved so that no kinematic singularities 
are introduced in the process, i.e., all gj remain regular. 
In analogy to Eq. (7) for the two-photon current example, 
one must choose the pj such that they carry no kinematic 
prefactors. In practice this is not always possible: there 
are equations where all pj come with kinematic prefac¬ 
tors and one must divide by them, thereby introducing 
kinematic singularities. Since the pj are regular, some of 
them must vanish in these kinematic limits. Therefore, 
the division should be done such that only the minimal 
number of pj picks up kinematical zeros. 

After reinserting Eqs. (119) into the general expres¬ 
sion (116), the resulting amplitude will take the form 

41 95 

= Y,fi rpr + E 9= (120) 

which is the analogue of the two-photon current (13). 
The first term is the transverse part of the amplitude, 
with transverse tensors that have mass dimension 

4, 6 , 8 ..., and dressing functions fi that become constant 

in any kinematic limit. The second term constitutes the 
gauge part, which is neither longitudinal nor transverse. 
Here we have again added the pj , which we eliminated in 
the first place; consequently, the gauge part must vanish 
if the amplitude is gauge invariant. In turn, if it does 
not vanish gauge invariance must be violated — either 
by a calculation that respects gauge invariance but is 
incomplete, or by an approach where gauge invariance is 
simply not built in. 

The fact that the Tj remain with mass dimension 
0 , 2 ,4,... is also the reason why violating gauge invari¬ 
ance can have severe consequences in practice. With an¬ 
other transverse projection of Eq. (120) everything col¬ 
lapses into the transverse part, in the same manner as 
in Eq. (14). If the dressing functions pj are nonzero, 
they will introduce artificial singularities with momen¬ 
tum powers -4,-2, etc. into the fi. In any case, the 


tP = P^ ■PjS'^-'-p^ Pi, 

— ptivaP a 8 
^ij • ^ r'i 


( 121 ) 


as the building blocks for the construction of such ten¬ 
sors [33]. tp is transverse with respect to p/ and pJ, and 
ep is transverse to both momenta. The only two tensor 
structures with mass dimension four are then 


^p8- = tptp and = ( 122 ) 

They have a simple physical interpretation: tp is the 
leading tensor of a scalar two-photon current with photon 
momenta pi and P 2 (which we now count as incoming), 
and eP is that of a pseudoscalar two-photon current (e.g. 
for the process tt — 77 ). Hence, if the LbL amplitude 
has scalar poles, they will appear in the form factor of 
ipi whereas the form factor of 1/2 inherits the pion pole. 

Next, we employ these tensors as permutation-group 
seeds in analogy to the derivation of Table IV. We take 
tpi and 'ip 2 as seed elements and derive the multiplets 
according to Eqs. (43-46). It turns out that each of them 
generates a singlet S and a doublet of type Vi ; the other 
multiplets vanish. The only singlets of dimension n = 4 
are therefore the tensors 

^ Irl ) — *12 *34 + *23 *14 + *31 *24 ! , . 

C7,/, 1 c-P'^ I c'^P _L c-PP ' 

O \PP2) — Si2 £34 + ^23 ^14 + £31 ^24 J 

and their corresponding (fully symmetric) form factors 
should be expected to be the dominant ones. Here and 
in the following we denote the multiplets for the Lorentz 
tensors with primes to distinguish them from the mo¬ 
mentum multiplets. 


We note that in the context of the LbL amplitude not even the 
constituent-quark loop is truly gauge invariant [21], Instead, the 
gauge part is a constant, + S’^P 5^°' + Spp 5‘"^)/(247r^), 

and drops out if the identity = —p\dMP'^P^Idp\ is 

employed as it is typically done in g — 2 calculations. 
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n 

Seed element 

# 

Multiplets 

n = 4 

n — 6 

n = 8 

n = 10 

n = 12 

4 

.pa- 
^12 ^34 

3 

5, Vi 

1 

1 

1 




pu pa 
^12 ^34 

3 

S, Vi 

1 

1 

1 



6 

aXa.au pX^.^a 

C22 ^3 ^44 

12 

S, !Di, 2^2, ^ 


1 

3 

5 

3 


,pX ,Xa 
^12 ^33 ^44 

6 

5 , Ui, T0 


1 

2 

3 



.pu ,pX ,Xa 
^12 ^31 ^24 

7 

s, r+, Tf 


1 

1 

3 

2 


p\ .Act 

=12 =31 '24 

7 

^2, 1 Ti ‘ T 2 



2 

5 


8 

,pu .pa .ap .0a 

^12 ^31 ^12 ^24 

3 

5, Ui, r+ 



1 

2 



Total 

41 


2 

5 

11 

18 

5 


TABLE V. 41-dimensional tensor basis for the transverse part of the photon four-point function, n denotes the mass dimension 
of the seed elements and # the number of the resulting singlets with mass dimension n. For the grayed terms we keep only the 
lowest-dimensional singlets with n = 10. 


To obtain the remaining basis elements, we define 


A 

• c. . 


(124) 


These are the lowest-dimensional Lorentz tensors that 
are transverse without introducing kinematic singulari¬ 
ties. is transverse to the momentum and is 
transverse in all Lorentz indices. Both are antisymmet¬ 
ric in a and /3.^^ The quantities in Eq. (121) are their 
momentum contractions: 


,/ii/ _ .fiav a 
‘'ij ‘'i Pj J 


= £r''p^ 


-ij '-i i'j ) 

and by contracting once more we can define 

lU ct t-f' U 

hjk ■= Pk=Pi- Pj Pk-Pi- Pk Pj , 
^ijk ■— ^ij Pk - ^ Pi Pj Pk ■ 


(126) 


(127) 


At dimension n = 6 one can find many possible Lorentz 
tensors by taking suitable combinations of these quanti¬ 
ties and their momentum contractions; however, only few 
of them are linearly independent. In particular, it turns 
out that the seed elements 


/ Ui/po- _ pAa ,C(1/ pA^ ,/ 3 cr 
rs — k-i (-22 ^3 ''44 J 

jM-vpcr _ ,p.v ,pA ,Act 

‘C4 “ ''12 ''33 '^44 : 

ipi^po- _ ,k.v ,pA ,Act 

Yb ~ ''12 ''31 ''24 : 

/ ui/po- _ ui/ pA , Act 
Y6 ~ ^12 ^31 ''24 


(128) 


Note that the electromagnetic field-strength tensor and its dual 
can be expressed in terms of these quantities: 

( 125 ) 


together with one element at n = 8, 

(129) 

are sufficient to generate a complete tensor basis with 41 
elements. The multiplets that they produce are collected 
in Table V. We will discuss them in a moment, but let 
us first resolve the remaining issue. 

Ultimately we would like to cast the 41 tensor elements 
into permutation-group singlets, so the question is: how 
can one construct singlets from, for example, a doublet? 
According to Eq. (50), the only possibility is to take dot 
products with other doublets. They will be made from 
the momentum multiplets iSq, T^o and To that we defined 
in Eqs. (84-86): 5o has mass dimension two whereas Vq 
and To have dimension four. The second row in Table III 
collects all possible doublets at dimension four and six. 
Combined with Dq, these are 


Vq , 'Do*'Do , To * To , 'Do* {To* To) (130) 


apart from further trivial multiplications with singlets. 
Now take for example the doublet D' — Vi{ipi), which 
is obtained from the tensor structure '0i. It has two in¬ 
dependent components, and by dotting it into two inde¬ 
pendent doublets from the list above we can generate two 
singlets. Restricting ourselves to the lowest-dimensional 
possible combinations, these are 


S[=Vo-V', 

S 2 = {a'Do *Do + PTo * To) • D', 


(131) 


where a, /3 are constants. Hence, in the process of con¬ 
structing singlets we have raised the dimension by two 
and four, respectively: a doublet with dimension n gen¬ 
erates a singlet at n -|- 2 and another singlet at n -I- 4. 
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To construct singlets from a triplet, we have to dot 
it into three independent triplets. The available options 
from Table III are 

To , To V ; To V To (132) 

which have dimension two, four and four, respectively. 
Therefore, a triplet of dimension n generates a singlet at 
n + 2 and two singlets at n + 4. Similarly, the fourth 
row in Table III shows that an antitriplet of dimension 
n produces a singlet at n + 4 and two singlets at n + 6, 
and the last row entails that an antisinglet of dimension 
n leads to a singlet at n + 6. 

In this regard. Table V should be read as follows. We 
start from the seven independent seed elements '0i ... '07 
defined above, with dimensions n = 4,6, 8. With their 
help we can generate 41 linearly independent tensors be¬ 
cause each seed generates a number of multiplets. Com¬ 
bining them with the momentum multiplets in the way 
described above, we generate further singlets whose di¬ 
mension has raised: the singlets have dimension n = 
4,6,8,10,12. 

Working out all multiplets by hand would be rather 
tedious, but it is easy to implement in Mathematica. We 
start from a large number of seed elements (essentially all 
conceivable tensor structures at n = 4,6,8 etc.) and let 
Mathematica generate the multiplets. We then add them 
up until all singlets at n = 4 are found, proceed to n = 6, 
etc. At each step we check for linear independence, i.e., 
whether the condition 

N 

if X! -C* q,k) = 0 => fi = 0 (133) 

i=l 

is still satisfied, until = 41 is reached (or N = 136 
in the case of Table IV). In that way one also confirms 
directly that there cannot be more than 41 linearly inde¬ 
pendent transverse elements (or more than 136 elements 
in general). 

Of course it is possible to construct many equivalent 
bases by this procedure, but they share some common 
features. First, the maximum number of singlets for a 
given mass dimension n (the last row in Table V), ordered 
with increasing n, is fixed: we found at most two singlets 
with n = 4, five singlets with n = 6, etc. Second, we 
never found fewer than five singlets with n = 12. By 
contrast, the same procedure applied to Table IV would 
produce singlets with n = 10 at most (they originate from 
the antitriplet and antisinglet seeds with n = 4 whose 
mass dimension is raised by 6.) The appearance of n = 12 
singlets in the transverse basis therefore suggests that not 
all of them are related in a simple way to the type-I basis 
without kinematic prefactors, as it was the case in the 
two-photon current example (8), and that divisions must 
have been necessary in the solution process of Eqs. (II9). 
Barring oversights, we are therefore led to believe that 
Table V can indeed serve as a minimal basis for the LbL 
amplitude. 


There is, however, a remaining problem. The construc¬ 
tion of singlets with lowest mass dimension from a given 
multiplet is not unique, as one can infer from the param¬ 
eters a, j3 in Eq. (I3I). There are two momentum dou¬ 
blets at n = 4 and both of them are equally suitable for 
constructing a singlet. Choosing one over the other can 
result, once again, in kinematic singularities. Similarly, 
there are three antitriplets at n = 6 and two antisinglets 
at n = 6. In fact, only the triplet case is unique because 
there are three momentum triplets up to n = 4 (those in 
Eq. (132)). One might conclude that it is simply impos¬ 
sible to construct a 41-dimensional transverse basis made 
of singlets, thus effectively leading to a redundant basis. 
On the other hand, one can argue that the solution of 
the system of equations (119), which we circumvented so 
far, should be unambiguous and determine these coeffi¬ 
cients in the process. (We mean ‘unambiguous’ in the 
sense that linear combinations of singlets with the same 
mass dimension are still allowed.) Ultimately it might 
turn out to be unavoidable to solve Eqs. (II9) directly, 
because even with a 41-dimensional transverse tensor ba¬ 
sis at hand one still needs to construct a gauge part that 
is kinematically safe and consistent with it. 

We were recently made aware^^ of a similar attempt 
in constructing a transverse basis for the LbL ampli¬ 
tude [40]. Expressed in our language, the seed elements 
(Eq. (3.14) of Part III therein) have the form 
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^13 

4 . 0 iiy\ 

^42 ) 

.poiX .ajSX 
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The problem of minimality is not addressed, but in terms 
of counting mass dimensions these tensors are equivalent 
to the seeds in Table V: there are two seeds with n = 2, 
four with n = 6 and one with n = 8. After working 
out the permutations, also the distribution of singlets 
is the same: we found an ‘optimal’ arrangement where 
41 = 2 -f 5 -I- II -|- 18 -f 5, as in the last row in Table V. 

Finally, we should comment on the four-gluon vertex. 
In that case the construction of the transverse part is the 
same as in Table V, except for the final construction of 
the singlets. The reason is of course the additional color 
structure which also produces multiplets (see Sec. IV C), 
so there are more possibilities for constructing singlets 
which have lower mass dimensions. In addition, the four- 
gluon vertex has a nonzero gauge part, which makes a so¬ 
lution of Eqs. (II9) mandatory. On the other hand, the 
presence of the gauge part also simplifies the problem be¬ 
cause a type-I basis such as that in Table IV (or even the 


We thank Gilberto Colangelo for bringing this to our attention. 
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one in Eq. (110)) would suffice for its solution. The con¬ 
traction with transverse gluon propagators (in Landau 
gauge) will remove the purely longitudinal elements, so 
that only the ‘survivors’ remain. This is exemplified by 
a calculation of the three-gluon vertex, including its full 
tensor structure and full kinematics, with a type-I basis 
only [7]. To reiterate, the problem of transversality and 
analyticity is only truly relevant when gauge-invariant 
amplitudes are considered. 

VII. SUMMARY AND CONCLUSIONS 

In this work we discussed in some detail the applica¬ 
tion of the permutation group S 4 to four-point functions 
of particles with four external gauge bosons. In particu¬ 
lar, we explored the case of light-by-light scattering that 
has important applications in the calculation of hadronic 
contributions of the anomalous magnetic moment of the 
muon. The main problem was to identify the transverse, 
gauge invariant components and represent these in terms 
of tensor structures that are free of kinematical singu¬ 
larities. To this end we made judicious use of the per¬ 
mutation group. We introduced an efficient notation for 
constructing the S 4 multiplets, and we applied it to or¬ 


ganize the kinematic phase space and to construct ap¬ 
propriate tensor bases. The multiplet analysis provides 
us with a straightforward way to find a complete basis of 
136 linearly independent elements. Implementing gauge 
invariance, we constructed a transverse basis of 41 ele¬ 
ments that has the required analyticity properties. The 
remaining problem is to cast this basis into one made of 
permutation-group singlets and combine it with a consis¬ 
tent gauge part. This will have to be addressed in future 
work. Nevertheless, we do hope that our general frame¬ 
work serves to establish common grounds for communica¬ 
tion among the several groups studying the light-by-light 
scattering amplitude in different approaches. 
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